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^ . In this paper we construct generalizations to spheres of the well known Levi-Civita, 

H ! Kustaanheimo-Steifel and Hurwitz regularizing transformations in Euclidean spaces of 

dimensions 2, 3 and 5. The corresponding classical and quantum mechanical analogues 

of the Kepler-Coulomb problem on these spheres are discussed. 
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1 Introduction 

It is well known that the problem of a body moving under the influence of a central force 
field with potential V{r) = —fj^/r has a singularity at the origin. We refer to this as the 
Kepler problem. This problem is usually posed in 3 dimensions, but since the motion is 
always constrained to a plane perpendicular to the constant angular momentum vector 
we can reduce it to 2 dimensions with Newtonian equations of motion and energy integral 
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where r"^ = r ■ r, r"^^ = c and r = {x,y) = (r cos 9, r sin 9). As is well known |jT|, |^, 
in two dimensions the Levi-Civita transformation effectively removes the singularity and 
rewrites this problem in terms of the classical harmonic oscillator. In this process the 
original problem has been regularized. To achieve the regularization, instead of t we use 

the variable s defined by 

dt did 

s 



dt r ds 



With x' 



^ etc., the original equations (^ are 



r' , a 
—r' + -r 



0, ^r' ■r'-^ = h. 
2r r 



Instead of using the variables {x, y) it is convenient to make the transformation 
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or 



L(u)u. 
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IFtoiu. the explicit form of these relations it follows that r' = 2L{u)u'. The equations of 
motion are equivalent to 



u 



u ■ u 



-u 



u, — 
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u ■ u 

Consequently we have the regularized equation of motion 

h 



u ■ u u ■ u. 
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This is essentially the equation for the harmonic oscillator if /i < 0. The solution ui = 
acos(cc;s), U2 = l3sm{us), up' = —h/2 is equivalent to elliptical motion. 

The relationship between the harmonic oscillator and the corresponding Kepler prob- 
lem can also be easily seen from the point of view of Hamilton- Jacobi theory. Indeed the 
Hamiltonian can be written in the two equivalent forms 



/^ 



H=-ipl+Pl) + ^^—-2 
2 ^ Vx^ + y' 
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(Mf + Mi)t^-+^-+'^]- 



(6) 



If we now write down the corresponding Hamilton- Jacobi equation via the substitutions 



Pni 



du,S 



j^i 



Pv 



(Ju2^ 



O,, 



we obtain 

Si + Si + 8/x - 8E{ul + ul) = 0. (7) 

This is just the Hamilton- Jacobi equation for a mechanical system with Hamiltonian 

H' = pI+pI-8E{uI + uI) 

and energy — 8/i. (This is the pseudo- Coulomb problem, see 0, |^. Reference Q also 
obtains (^ as an application of Stackel transform theory). 

This transformation also achieves a regularization of the corresponding quantum me- 
chanical problem, which we call the quantum Coulomb problem. Indeed, the Schrodinger 
equation in the presence of the potential V{r) = —fj^/r in two dimensions has the form 

In the coordinates (mi,U2), (||) becomes 0] 

(dl + dlj<^ + {8/i + SEiul + ul)}<^ = 0. (9) 

Here, (^ has all the appearances of the Schrodinger equation in a oscillator potential 
V{ui,U2) = —4:E{u1 + ul) and energy S = 4/i. Note that for scattering state E > 
we have the repulsive oscillator potential and for E = the free motion. For E < we 
get the attractive oscillator potential and the corresponding bound state energy spectrum 
can be easily computed from this reformulation of the Coulomb problem, although the 
weight function for the inner product is no longer the same 0, |^, ^ |^ . (Indeed, the Virial 
Theorem states that for the Coulomb problem the change in weight function does not 
alter the bound state spectrum, ||^). The wave functions have the form $ = (pi{ui){p2{u2) 
where the functions ipx satisfy 

{dl + KX + 8Eul)ipx = 0, A = 1, 2, ki + K2 = 8fi. 

The bound state eigenvalues are quantised according to 

Kx = 2V^2E{2nx + l), A = 1,2, (10) 

where ni,n2 are integers. Taking into account 0, $(— Ui, — ^2) = (— l)"'^^"^$(ni, U2) and 
using that \E'(x) is even in variable u: \E'[a;(n)] = \E'[2;(— m)], [because two points (— Mi, — M2) 



and {ui,U2) in -u-space map to the same point in the plane {x,y)] we find /^from ([T0|) the 
energy spectrum of the two dimensional Coulomb system [H, 0] 



/U^ , . rii + n2 



^"'-WW ^=^ = 0-l-2 



\2' 



It is well known that the regularizing transformations @) that we have discussed for 
the Kepler and Coulomb problems in two dimensional Euclidean spaces are also possible 
in the case of three (Kustaanheimo-Stiefel transformation for mapping R4— s^Rs) p|, |l^, 
[TT], |T^ and five (Hurwitz transformation for mapping Rg— s^Rs) [0, |I^, |I^, 0, |T^, [T^, 
p!9| dimensions. The only difference in these cases is that additional constraints are 



required. These transformations have been employed to solve many problems in classical 
and quantum mechanics (see [|I^] and references therein). 



As in flat space, the study of the Kepler- Coulomb system in constant curvature spaces 
has a long history. It was first introduced in quantum mechanics by Schrodinger |]2D|, who 
used the factorization method to solve the Schrodinger equation and to find the energy 
spectrum for the harmonic potential as an analog of the Kepler-Coulomb potential on 
the three-dimensional sphere. Later, two- and three-dimensional Coulomb and oscillator 



systems were investigated by many authors in [^ |^, |2^, |2J, ^, ^ ^7\, ^ |^, |3^, |3T| . 

However, in spite of these achievements the question of finding all transformations that 
both generalize the Levi-Civita, Kustaanheimo-Steifel and Hurwitz transformations for 
spaces with constant curvature and preserve the Kepler-Coulomb and oscillator duality 
has been open till now. The answer to this question is a main aim of our paper. 

The paper is organized as follows. In §2 we present the transformations that generalize 
the flat space Levi-Civita transformation and correspond to the map S2c~^5'2 /,from com- 
plex into real two-dimensional spheres. We show also that this transformation establishes 
the correspondence between Kepler- Coulomb and oscillator systems in classical and quan- 
tum mechanics. In §3, in analogy with §2, we construct the Kustaanheimo-Steifel and 
Hurwitz transformation and show Kepler- Coulomb and oscillator duality for mappings 
S4C— ^Ss and Ssc— ^-Ss, respectively. §4 is devoted to a summary and discussion of our 
findings. In the Appendix we give some formulas determining the connections between 
Laplace-Beltrami operators and the volume elements in different spaces. 

2 The transformation on the 2-sphere 

The potential, which is the analogue of the Coulomb potential in quantum mechanics and 
the gravitational potential for the Kepler problem, is taken to be EO, |2TI 



V=-^^^^, (11) 

where {si, 82,83) are the Cartesian coordinates in the ambient Euclidean space and R is 
the radius of the sphere 

s ■ s = sl + sl + sl = R"^. 

(Note that V = — ^ cot a where a is the arc length distance from s to the north pole of 
the sphere. Furthermore, the leading term in the Laurent series expansion in a about the 
north pole is — -^)- 

This problem is easily transformed into a much simpler one via the transform 

9 2 



Si = i\/ul + U2 + ul 



2u. 



S2 = i\^ui + ui + u% , (12) 

«3 



,2 I „,2 I „,2 |„, I "^1 + "^2 



■S3 = V"l + "2 + M3- M3 + 



2m, 



or in matrix form 



S3 



Ui + U2 



Un 



2m, 



iui —iu2 
iu2 iui 
Ui U2 2m3 



Ml 
U2 

U3 



(13) 



The advantage of this transform is the Euler identity 



2\2 



+ S2 + S3 = (Ml + M2 + %) 



(14) 



^from which we see that the point u = {ui,U2,u^) hes on the complex "sphere" '^20'- 
ul + u^ + u^ = D"^ with the real radius D ii s = {si,S2, S3) lies on the real sphere S2 with 
radius R, and R = D"^. 

In the general case the two-dimensional complex sphere 820 Kiay be parametrized by 
four real variables (the constraint u{ + u\ + u\ = D^ includes two equations for real and 
imaginary parts). The requirement of reality of the Cartesian variables Sj leads to two 
more equations and the formula (0) corresponds to the mapping /^from a two-dimensional 
submanifold (or surface) in the complex sphere S2C (four dimensional real space) to the 
sphere S2. To verify we introduce ordinary spherical coordinates on S2: 



si = R sin X cos 99, S2 = R sin x sin y?, S3 = R cos x- 
iFroTsi transformation (0) we have 

R~ 2\D M3 



(15) 



(16) 



Putting S3 = R cos X in formula 
on the complex sphere 820 are 



we get U3 = De^^ and then the corresponding points 



Ml 



dVT 



o2ix 



COS ■ 



^ 



U2 



dVi 



o2ix 



sm ■ 



V 



M3 = De^^ 



(17) 



where < X ^ '^^ 0<V9< An. Note that the transformation (|T2D is not one to one; two 
points (—Ml, — M2, M3) and (mi, M2, M3) on the sphere in M-space correspond to one point on 
the sphere in s-space. Thus, when the variables (mi,M2,M3) cover the sphere in M-space, 
the variables Sj cover the sphere in s-space twice. 



Let us now introduce nonhomogeneous coordinates according to |^ 



R- 



Ui 



D- 



Ui 



D' = R, 



S3 



U3 



1,2. 



Then formula (|T^ transforms to 



Sl 



i{ul - ul) 



2(1 + 



2D2 ) 



S2 



IU\ ■ M2 



:i + 



2D2 > 



(19) 



In the contraction limit D — i> 00 we obtain 

Sl = i- 



.u\-ul 



S2 = iUi ■ M2, 



(20) 



which coincides with the flat space Levi-Civita transformation (D up to the additional 
mapping Ui -^ e~^^\/2ui. 

The relationship between the infinitesimal distances is 



ds ■ ds = {ul + ul + ul) 



{u-duY ful + ul' 



du ■ du 



ui 



Un 



3(n ■ duY 



Thus, when restricted to the sphere, the infinitesimal distances are related by 



ds ■ ds (u^ + 'uVX , 

— =; — = — 5 — du ■ du, 

R \ ui 



and we see that as in flat space the transformation (|T2D is conformal. 



(21) 



:22i 



2.1 Classical motion 

Just as in the case of Euclidean space, the classical equations of motion under the influence 
of a Coulomb potential can be simplified. The classical equations are 



s = -(s ■ s)s- W, 



(23) 



where the first term on the right hand side is the centripetal force term, corresponding to 
the constraint of the motion to the sphere, and the potential satisfies 



s-W = 0. 



(24) 



Here, s = -^s. (In studying (|23D and (|2J) we initially regard the coordinates s as un- 
constrained and then restrict our attention to solutions on the sphere). In the case of 
potential (111]) these equations become 



d . . . ^ /i SjS^ 

-—t:Sj = —Sjls ■ s) — — — 

dt' ' '^ ' R{sl + sl 



J = 1,2 



d"^ yU 1 

TI^SS = -S3 S • S) + -— -J 

dV Ri^sl + sD'i 



subject to the constraints 



s- s = R^ 

and its differential consequences 

iFiova. the equations of motion we immediately deduce the energy integral 

-s- s + V = E. 
2 

We choose a new variable r such that 

dr 1 m| 

m^^' ul + ul' 



(25) 



(26) 
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In terms of the variables r and Ui, the equations of motion can now be written in the 
form 



2/i 



2D^ 



2/i 



u[y + {u',y + {u',Y - 2D\E + -^) + ^(E - -^ 



D2 



u% 



D^ 



u'l + 2{E + ^)u, = 0, u'i + 2{E+'-^)u^ 



Z/i 



2D^ 



Z/i 



"S + 2(S+P"3-^(B-^) 






(27) 





(28) 


0, 


(29) 



subject to the constraint u ■ u = D^ and its differential consequences u ■ u' = Q, u ■ u" + 
It' ■ w' = 0, where u[ = ^. These equations are equivalent to the equations of motion we 
would obtain by choosing the Hamiltonian 



H = \{vl+vl,+vl)-{E+'J^){ul + ul + u^,) + ^^{E-'J^^, 



D^ 



D^ 



(30) 



regarding the variables Ui as independent and using the variable r as time. In fact, to 
solve the classical mechanical problem from the point of view of the Hamilton- Jacobi 
equation, we use the relation 



\{pi,+pi^ + p. 



/i S3 



R'Isl + sl 



E 



Un 



Ui + U2 



1 ' 2 , 2 , 2 N !, ^^ , t:^\ . D\t:^ ^/^ 



K<+<+<)-(^^ + ^) + %(^-^) 



2^2V^«i ^"2 ^u„ V ^2 ' ' • ui^ D^' 



(31) 



together with the substitutions p^. = ^ and Ps^ = ^ to obtain the Hamilton- Jacobi 
equations 



dsi. 



+(^]\(^y 2/i 



S3 



ds2j \ds3J 



R 'Isl + sl 



2E = 0, 



dS 



'dS 
duo 



+ 



' dS 



^9,m , 2D^ , iu 

2D^(^ + E) + ^{E ^ 



'D2 



Un 



D2 



0. 



(32) 



(33) 



\dUiJ \UU2J \L/"3/ ^ "3 

This last equation can be solved by separation of variables in the spherical coordinates 
on the complex sphere S2C (|T7|) . 



2.2 Quantum motion 

If we write the Schrodinger equation on the sphere for the Coulomb potential (pJ]) 



1A.VI.+ (E + ii_iL 
2 \ R.^sl 



^ = 0, 



(34) 



and use the transformation (12), we obtain [see formula (|145|) ] 



1^ , / u'^D^ul + uW , , 



(35) 



where 



E = 2i/i, 



UJ 



2(.-^) 



(36) 



Thus we see that the Coulomb problem on the real sphere S2 is equivalent to the cor- 
responding quantum mechanical problem on the complex sphere S2C with the oscillator 
potential (Higgs oscillator [^, |26|, ^) and energy 2i^, but with an altered inner product 
(see the Appendix). 



Let us consider the Schrodinger equation (35). Using the complex spherical coordinates 
(|T^) we obtain 



Id dtp 1 dip 

smx^r + ^^—TT 

sm X ox ox sm x o^ 



+ uj^D^ - iSD' 



3«X 



smx 



ij = o. 



(37) 



S = 2k. 



UJ 



2(^-7^)- 



(3^ 



To solve this equation we first complexify the Coulomb coupling constant /i by setting 
k = ifi in the formulas for S and u 

k_ 

Further, we analytically continue the variable x iiito the complex domain G: < Rex < tt 
and < Imx < cxd, (see Fig.l) and pass from the variable x to i), defined by 

e^x = cos^. (39) 



Im X 



G 



G 



TT 



Rex 



Figure 1: Domain G ={ < Re x ^ tt; < Im x < 00} on the complex plane of x- 



For real 9 this substitution is possible if Rex = or Rex = tt and Imx € (0, oo), which 
corresponds to the motion on the upper (0 < ^ < |) or lower (| < "(9 < vr) hemispheres of 
the real sphere. In any case conditions (|39|) and (|38|) translate the oscillator problem from 
the complex to the real sphere with spherical coordinates {'d,(p/2). In these coordinates 
we can rewrite ( |37D in the form 

— -— sin^-^ + -^TT^ + i2SD^ + uj^D^) — U = 0. 40 

sm i!) av av sm d OLp^ \ cos^ v J 

Using the separation of variables ansatz 

i,{d,^) = R{df^, m = 0,±l,±2,... (41) 



we obtain 

sm V dv dv y cos^ v sm v J 

The corresponding solution regular at the points -(9 = 0, 7r/2 takes the form 0] 

Rnrmi-d) = Cn,.m{i^) (sm^y""^ {cos^y^KFi{-nr,nr + u + \m\ + 1; |m| + l;sin2^) 

= Cr^^rniiy) }''^]'^!^L (sin^)!-' (cos ^)-+^ /^(H'^) (cos 2^) (43) 

with energy spectrum given by expression 

S = ^[(n + l)(n + 2) + {2u - l){n + 1)], v = (^'^' + ^)' (44) 

where C„^m(i/) is the normalization constant, P'^"'^-'(x) is a Jacobi polynomial, rir = 
0, 1, 2... is the "radial" and n = 2nr + \m\ is the principal quantum number. 

To compute the normalization constant Cn^mi^) for the reduced system we require 
that the wave function (^) satisfy the normalization condition (see the Appendix): 

- — / ^n^mrn^m^^ dv{u) = D' / Rn^^Rl^ siu X^X = 1 (45) 

2 JS2C ""3 Jo 

where the symbol "o" means the complex conjugate together with the inversion x -^ —Xi 
i.e. tp^ix^ (p) = ip*{—Xy '{')■ [We choose the scalar product as ip'^ip because for real tu^ and 
£ the function V'*(X) ^) ^-Iso belongs to the solution space of ([371).] 

Consider now the integral over contour G in the complex plane of variable x (see Fig.l) 

//■TT /'7r+200 

^n,mK.m sin XC^X = j RurmRl^m^^T^XdX + J RurmRl^m^^^XdX 

rioo rO 

+ / K^^i?^^^ sin xrfx + /. K.mK.m sin xc?X- (46) 



Using the facts that the integrand vanishes as e^*'^^ and that -R„^m(x) is regular in the 
domain G (see Fig.l), then according to the Cauchy theorem we have 



/•TT rioO /'TT + ZOO 



1 _ e^i^i-^^^) 



^nrm^rirm ^^^ X'^X- 



Making the substitution ( [T7| ) in the right integral of eq. (^ , we find 



^rirmJ^nrm ^^^ X'^X 



I _ g2i7r(i.+ l) 



[i^^^J^sin^tan^^rf^. 



(47) 



(4J 



Using the following formulas for integration of the two Jacobi polynomials ||33| 

' {I - xY{l + xf~^[P^''^^\x)fdx 
1 

we find 



2"+/3+ir(n + a + l)r(n + /? + !) 
(2n + a + /5 + l)n\T{n + « + /? + !) 

2"+/3r(n + a + l)r(n + /9 + 1) 
(/5)n!r(n + a + /3+l) 



-z/(z/ + 2nr + \m\ + 1) (n,, + |m|)!r(|r7i| + n^ + z/ + 1) 
(|m|)! ^ D4[l-e2*'^('^+5)](2n^ + |m| + l)(n^)!r(n^ + i^+ 1) 



(49) 



The wave function Tpiy'd, (p) = Tpnrm.{'&, ^) is then given by eqs. (0), ( |i3| ) and 

Now we can construct the Coulomb wave functions and eigenvalue spectrum. From 
transformation 



^nrm{^, yp + 27r) = e'^'">„,^(^, ip) 



(50) 



and requirement of 2n periodicity for the wave functions (0) we see that only even 
azimuthal angular momentum states of the oscillator correspond to the reduced system. 
Then, introducing new angular and principal quantum numbers M and A^ by the condition 

n = 2nr + \m\=2nr + 2\M\ = 2N, N = 0,1,2,..., |M| = 0, 1, 2, ..A^, (51) 

comparing (|38D with expression (^4]) for the oscillator energy spectrum, and putting k = 
ip,, we find the energy spectrum for reduced systems 



E 



N 



N{N+1) 
2i?2 



/i- 



2(iV+i)2- 



(52) 



This formula coincides with that obtained from other methods in works pT] , p6| , ^% 
Transforming {} back to the variable x by (11), we see that (|44D and (|3^) imply 



iy = icT-{N + -), 



a 






1 ■ 
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Using 



r(l/2+|M|+za) 
T{l/2-\M\+ia) 



^ |r(l/2 + ^a)|2 



.niMi 



TT 



cosha7r|r(l/2+ |M| +ia)p 



we easily get from (|4l|), (|43|) and ( ^9]) the eigenfunction of Schrodinger equation 



^7vm(x,<^) = C;vAf(^)e-*^(^"l^-'l-^'^)(sinx) 



M| 



X 2i^i(-A^ + \M\, \M\ + «cr + i; 2|M| + 1; 1 - e^'^' "" 



iMip 



2tx 



(53) 



(54) 



where now 

CArA/(o") 



2\M\ 



[(iV+i)2 + a2](iV+|M|)! ^ 



i?(2|M|)!\ vr(iV + i)(iV-|M|)! 



e-|r(|M| + l/2 + za)| 



(55) 



By direct calculation it may be shown that the Coulomb wave function 
normalization condition 



satisfies the 



/■TT r-ln 

R^ / sinxc^X / d^-<^NM'^*N'M' 
Jo Jo 



^NN'^AIM'- 



Thus, by reduction from the two-dimensional quantum oscillator on the complex sphere 
we have constructed the wave function and energy spectrum for the Coulomb problem 
on the two-dimensional real sphere 5*2 . Formula (^^ for Coulomb wave functions on the 
two-dimensional sphere is new. 

Now let us consider the flat space contraction. In the contraction limit i? — *■ oo the 
energy spectrum for finite N goes to the discrete energy spectrum of the two-dimensional 
hydrogen atom |^, || 



lim En(R) 



f^^ 



2(iV+i)2' 



A^ = 0,1, 



~ V ~ ^. where r is the radius- vector in the two- 



X 



In the limit R —>■ oo, putting tan^ 

dimensional tangent plane and using the asymptotic formulas p4 | 



i?' 



lim 2Fii-N + |M|, \M\ +ta + -; 2\M\ + 1; 1 - e^'^) 

R^oo 2 

X^O 



iF^{-N+\M\,2\M\ + l; 



2fir 



lim \r{x + iy)\e'2'^ \y\'2~ 

|j/|->oo 



^2tt, 






(56) 



11 



we obtain the well known Coulomb wave function with correct normalization factor 



lim ^atmIx^V^) 



(Ar+|M|)! / 2/ir Y^^' e ^ 
^^^ {N+\fl^\ N-\M\)\ \nT\J JmW- 



X iFi(-iV+|M|,2|M| + l;-^)-=. (57) 

In the case for large R and A^ such that A^ ~ kR, (where k is constant) we obtain the 
formula for continuous spectrum: E = k^/2. Now taking into account that o" ~ f and 
using the asymptotic relation (^), we have 



/k (2kr)^^^^ 



-ikr 



ill 1 P^^'f 

X iFi(|M| + f + -;2|M| + l;2^fcr)^, (58) 

fc Z v2vr 

which coinsides with the formula for the two-dimensional Coulomb scattering wave func- 



tion in polar coordinates [^ 



3 The three and five dimensional Kepler - Coulomb 
problems 

In complete analogy with the three- and five-dimensional Euclidean case the corresponding 
regularizing transformations exist for the Kepler and Coulomb problems in spheres of 
dimension 3 and 5. Indeed if we consider motion on the sphere of dimension n then the 
classical equations of motion in the presence of a potential are just (p3|), (|2^ ) again, where 
now 

s = (si,---,s„+i), (59) 

subject to the constraints 

s-s = R^ (60) 

and its differential consequences 
If we choose our potential to be 

these equations assume the form 

^^. = -^.---^^^7^, J = !,•••,-, (62) 

12 



^^ • • , At 
-— s„+i = -s„+is ■ s + — -y. (63) 

The energy integral again has the form (^). 

We are particularly interested in dimensions n = 3,5. We deal with each of these cases 
separately. 

3.1 Generalized KS transformation 

For n = 3 we choose the Uj coordinates in five dimensional space according to 



si = i^ul + ul + uj + ul + ul 



M1M3 + U2U4 
M5 



• / ~, 2~, 2~, ~, 2 ^2^^3 - U1U4 



S3 = iyul + ul + ul + ul + ul 



- uj + ul-ul-ul 
2uk 



5 



S4 = ^ul + ul + ul + ul + uhL + '''^''''^j'^'''A . (64) 



The basic identity is 



2, 2, 2, 2 /2| 2, 2, 2, 2\2 

Si + S2 + S3 + S4 = (Mi + ^2 + % + M4 + Mg) , 



and the basic relationship for the infinitesimal distances is 

dsl + dsl + dsl + dsl = (65) 

jii^l + ul + ul + ul) [dul + du\ + du\ + du\ + du^ 

+ {Uiidu^ — u^du^ + U2dui — Uidu2) }, 

where the constraint for mapping between the 3-sphere: I]i=i sf = -R^ and the complex 
4-sphere: Yl\=i uf = D^ is clearly 

u^duz — u^dui + U2dui — Uidu2 = 0. (66) 

In this section we will use the Eulerian spherical coordinates on the complex 4-sphere S4C 

Ml = DVl - e2^x cos f cos ^, M2 = Dy/l - e2^x COS f sin ^, 

M3 = D VI - e2»x sin f cos ^, M4 = /^Vl - e^^^sin f sin ^, (67) 

Ms = De^-^, 
where the ranges of the variables are given by 

< X < 71", < /5 < vr, < a < 27r, < 7 < 47r. 
The corresponding spherical coordinates on 5*3 are 

si = R sin X sin /3 cos a, $2 = R sin x sin /3 sin a, 
S3 = i?sinxcos/3, S4 = i?cosx- 
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3.1.1 Classical motion 

In analogy with our previous analysis we choose a new variable r according to 

dr 1 wi 



dt D^ ul + ul + ul + ul' 
In the u coordinates the equations of motion can be written as 

{u\f + {u',f + {u',f + {u',f + {u',f - 2D\E + ^) + ^-^{E - §) = 0, 

«;' + 2(E + -^)n, = 0, J = 1,2,3,4, (68) 

subject to the constraints 

fc=i fc=i 

5 

^{uku'l + {Uf.Y) = 0, M4M3 — M3M4 + M2M1 — M1M2 = 0. 
fc=l 

Note that equations (^) are compatible with these constraints. Here, the Kepler problem 
on the sphere in three dimensions is equivalent to choosing a Hamiltonian 

H = ^{pI,+pI,+pI,+pI,+pIJ (69) 

-(E + ^)iul + ul + ul + ul + ul) + ^{E - ^), 

regarding the variables Uj as independent and r as time. The only difference is that there 
is now the constraint 

UaPus - U^Pui + U2Pui - UiPu2 = 0. 

In terms of the Hamilton- Jacobi formulation we have the relation 

1(„2 ^2 ^2 ^ 2 ^_ /i _ ^4 _ p. ^ 



2 



«5 



uf + u| + M^ + uj 



1 



:(P«1+P«2+P«3+P«4+P«5 



L2D2 

D2 in. 1 

0. 
With the usual substitutions, the corresponding Hamilton- Jacobi equations are 



-(^ + ;^)(«i + 4 + ul + ul + ul) + -^{E - ^ 



S©^-(^-t7:rTlni-' 
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or 

and the constraint has become 
where operator C is 

C = U2' 



+ 



M5 



Z/i 



tfl 



-(E-^)-iE + ^^) 



d 



C-S = Q 



d d 

Ui— h U4- 



U3- 



d 



(71) 
(72) 
(73) 



' dui ^ du2 ' * du^ " du4 

Equation (^) can be solved by separation of variables in the spherical coordinates (|6 
on the complex sphere S^c- 

3.1.2 Quantum motion 

The associated quantum Kepler-Coulomb problem on the sphere corresponding to the 
potential (|6lD 



^a13)^+ \E + 



/i 



sa 



^-'4 + 4 + 4, 



^ = 



(74) 



translates directly to [see formula ( |152| )] 

2 " V 2 ul I 



with the constraint 
where C is given by (173|), 



and 



^ = 2^/^ - -]yi- 



£•$ = 0, 



\1/ = ^52$, 



uo^D^ = 2ED^ - 2ifi + 



4^2 



(75) 
(76) 

(77) 

(78) 



Here A^^-* and A^'^^ are Laplace-Beltrami operators on the spheres S3 and 84(7, respectively. 
Consider the Schrodinger equation (|75D in complex spherical coordinates (|67D. We 
have 



^-^— -- e '^ sm Xtt 
sm X OX OX 



^*x 



Lj'D'' - iSD' 



smx sm X 



$ = 



(79) 



where the operator L^ are defined in ( |150| ). We complexify the angle x to the domain 
G (see Fig.l) by the transformation (|39|) , such that {} G [0, |] and also complexify fi 
by setting k = ifi in expression for S and uj"^. Then equation (|7^) transforms to the 
Schrodinger equation for the oscillator problem on real sphere S4. 



15 



We make the ansatz 



where 



$(^,a,/3,7) = (sin^)-2 Z(^)P^^,^^(a,/3,7) 



^^„„..(«, A7) = e™^"C^,(/?)e^™^^ 



(80) 



is the Wigner function [^, satisfying the eigenvalue equation 

L^P^^,^^(a,/?,7) = £(£ + l)P^^,^^(a,/3,7), 
and normahzation condition 



Then the function Z{i)) satisfies 



d^^ 



+ 



2SD' + uo'D'' + 



9\ uo^D^ (2£+l 



f82l 



^3) 



,2 1 



z = o. 



4) 



4 J cos^ -(9 sin^ t? 

The corresponding solution regular at ^ = 0, 7r/2 and energy spectrum are given by 

ZnA^) = const (sin^)2^(cos'^)'^+5 2Fi(-n^, Ur + 2i + u + 2; 2i + 2; sin^^?), (85) 



S 



2D' 



:[(n + l)(n + 4) + (2z/-l)(n + 2)], 



(86) 



where u = iu'^D'^ + 4) ^) '"■ = 2nr + 2i = 0, 1, 2.. is the principal quantum number. The 
other quantum numbers are 

nr = 0,l,..n, 2£ = 0,l,...n, mi, m2 = —£,—£ + 1, ....,£— 1,£. 

Thus the wave function $(-(9, a,/?, 7) normalized under the condition (see Appendix) 



iD' 



27r Jsic 



dv{u) 






1 



^7) 



has the form 



/2£+ 1 

'^n.fmim2(^,a,/5,7) = C'n.K'^) U ~^^ ^"'■^(^) ^^i.^na («"^' 7) 



(88) 



with 



\i/+i 



i?^^^(^) = (sint9)^'(cos?9)^+5 2Fi(-n^, n^ + 2£ + z/ + 2; 2£ + 2; sin^t?), (89) 



Cnr(.{^) — —J 



[{-w){v + 2£ + 2nr + 2)](2£ + n,. + l)!r(2£ + z/ + n^ + 2) 



/}! \ (l-e2--)(£ + n, + l)[(2£+l)!]2(n,)!r(i/ + n, + l) 



(90) 



16 



We now construct the wave function and energy spectrum for the Schrodinger equation 
(|7^). The corresponding wave function \l/(s) connecting with $(m) by formula ( [77| ) is 
independent of the variable 7 and 27r periodic in a (the transformation a — > a + 27r is 
equivalent the inversion Ui — > — Mj, i = 1,2,3,4). The constraint ( [7BD in the spherical 
coordinate ( ]67| ) is equivalent to 



^7 



$n,£mim2 (^, ", ^, l) = "^2$n,£mim2 (^, ", /?, T) = 



and we have 7tt,2 = 0. From 27r periodicity we get that i and mi are integers. Then, 
upon introducing the principal quantum number N = {rir + £) + l = | + l and using the 
expression (|78D, we obtain the energy spectrum of the reduced system 



E 



N^ 



fX^ 



2/?2 



2m'- 



N=l,2. 



(91) 



where k = ifi. This spectrum coinsides with that obtained from other methods PO, El, 
Returning from i) to the variable x? observing that 



u = ia — N, 
and using the relations {mi = m) 

P^>,/?,7) = (-l)"^1 



a 



fxR 



4:71 



2£+l 



Yi-raiP, a) 



(92) 



we obtain the wave functions (with correct normalization) for the reduced system in the 
form 



^ 



Nim\ 



{x,P,a) = /De*^/^ ^nremoix,a,P,7) 



;-ir 



C„.^Acr) (sin a) e 



i -ia(N-i-ia) 



X 2Fi{-N + £+l,l + i + za; 2i + 2; 1 - e"''^) YimW, a), 



where 



C„e(a) 



2 ^ e 2 



{m + a^){N + i)\ 



T(l + i + ia)\. 



(93) 



(2£ + l)!^ 27riV(iV-£-l)! 
This solution is identical to that given for the Coulomb eigenfunction on S3 in papers 



23, 25]. Note that in [25| it already has been shown that the function (p3D contracts as 



R ^ 00 into the fiat space Coulomb wave function for discrete and continuous energy 
spectrum. 
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3.2 Generalized Hurwitz transformation 

The analogous problem in five dimensions can be realized via the variables 



Sl 



S2 



S3 



S4 



S5 



Sq 



J2^l] ~ (^1^5 + U2UG - u^ur - UiUs) , 

1 

{uiUq — U2U5 + UsUs — UiUj) 



Ug 

i 

Ug 



\k=l / 
f 9 N 

\k=\ J 

^9 \ 5 ■ 

y]"fc — (M1M7+M2M8+M3M5+M4M6), 

Kk=X I «9 
^9 \ I ^ 

V Mfc (wiMg - M2M7 - M3M6 + M4M5) , 

.^1 / «9 

^9 \l • 

V^ 21 '/2|2|2|2 2 2 2 2 

1 
f 9 

2_] Mfc I I "9 -r 77— Z^ "fc I , 



(94) 



^fc=l 



9fc=l 



which satisfy 



E-l= E 



M^ 



(95) 



The relation between the infinitesimal distances on the five dimensional sphere S5: 
Y^i=xs\ = R^ and eight dimensional complex sphere Sgc: Hi=iuf = -D^, {R = D"^ 
is 

1 A . o -1 



oE^^? = -r 



-Rj^"" ^ M5 



S "fc) E ^"^ + ^1+^2 + ^: 



9 L fc=i i=i 



(96) 



where 



tui = u^dui + Uzdu2 — U2du3 — uidui — u^du^ — u-jdu^ + ugduj + u^dus, 
UJ2 = u-^dui — Uidu2 — Uidu^ + U2du4 — Ujdu^ + ugduQ + u^duj — u^dug, 
^3 = U2dui — Uidu2 + Uiduz — u^du/i, + u^du^ — u^du^ + u^duj — u-jdus, 

and the constraint for mapping Ssc —^85 corresponds to 

coij = 0, "^ = 1, 2, 3. 

Following [|16| (see also |T8|) we can supplement the transformation (UM with the angles 



an 



Ph 



arctan 



2UiU2 

u\-u\ 



+ arctan 



2U3M4 



Mo — Ua 



E [0,27r) 



/y2 _j_ ^2\ 2 



(97) 
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iH 



arctan 



2uiU2 



— arctan 



2M3M4 
ul-u\ 



e [0,47r). 



The transformations ( p^ ) and (p7| ) correspond to 880-^88=85(8)83. If we now choose the 
spherical coordinates on 85 as 



Si + is2 = -R sin X sin -(9 cos |e* 
S3 + is4 = i? sin X sin ■& sin |e*' 



Q+7 



S5 = i? sin X cos "i?, 
Sq = Rcosx- 



then the corresponding (nonorthogonal) spherical coordinates on the eight dimensional 
complex sphere take the form [D"^ = R) 



Ml 
U2 
U3 
U3 
U5 



^ (3h OtH + IH 



ue 



Dvl — e^ cos — cos — cos ■ 

2 2 2 

n /I 2i^ ^ Ph . aH + lH 

JJV 1 — e"^*^ cos — cos — sm 

2 2 2 

n /I 2i^ ^ . PH an-lH 

Dy 1 — e"'*'<^ cos — sm — cos 

2 2 2 

n /I 2i^ ^ • -^^ • ^H-lH 

Uy 1 — e"'*'<^ cos — sm — sm 

2 2 2 

n /I 2i7 • ^/^ ^ ^^ a + 7 + «H + 7H 

iJ V i — e"'*'<; sm — cos — cos — cos 

2V 2 2 2 

. 13 . pH a--f-aH + lH 

+ sm — sm — cos 

2 2 2 

n /I 2i^ • ^/^ -^ /5h . a + 7 + a/f + 7/f 

Dyl — &^^^ sm — cos — cos — sm 

2V 2 2 2 

. (3 . Ph . a-j-an + lH 

— sm — sm — sm 

2 2 2 

n /I 2i7 • ^f- (^ f^H a--f + aH + lH 

DVl — e"'*^ sm — (sm — cos — cos 

2^ 2 2 2 

13 . pH a + 'j-an + lH. 

— cos — sm — cos ) 

2 2 2 ^ 

n /I 2i7 • ^^ • ^ ^^ • a-7 + aH + 7/f 

iJv 1 — e"'*^ sm — sm — cos — sm 

2^ 2 2 2 

P . Ph . a + '^ - aH + lH. 
+ cos — sm -— sm ) 

where x ^ [0,7r], ^ G [0,7r], a E [0,2tt], P G [0, tt] and 7 G [0,47r]. 



M7 



1*8 



(9^ 



3.2.1 Classical motion 

The Kepler- Coulomb potential on the five dimensional sphere 85 has the form 

Sfi 



--I 



sl + sl 



S3 + S4 + S5 



(99) 
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As before we can define a new coordinate r such that 

dr 1 Mg 



The corresponding equations of motion are given by 



EK)^-2(i? + ^)-^(i5-|^) = o, (100) 

e=i jy Ug jy 

ul + 2{E + ^) = 0, k = l,...,8 
„ , ill ^ 1D^ , in, , 

«9 + 2(i^ + ;^)«9-^(i?-;^)=0, 



subject to the constraints 

9 9 9 

Y^u\ = D\ Y.^,v!, = ^, ^(w,< + K)2) = 0, 

l.=\ l=\ £=1 

U^u'i + U3U2 — 'U2W3 — Uiu'^ — Ugu'^ — UjUq + Uqu'y + u^u'g = 0, 

U^u'i — U4U2 — Uiu'^ + U2U'^ — U-ju'^ + U^u'q + U^u'j — UQUg = 0, 

U2u'i — Uiu'2 + U4U'^ — U^u'^ + Uqu'^ — U^u'q + U^Uj — Uju'^ = 0. 

These equations of motion are equivalent to what we would obtain by choosing the 
Hamiltonian 

H=lJ:Pi-iE+^^)j:uj+^iE-^^), (101) 

regarding the variables Ui as independent and using r as time. The associated constraints 
are 

M4P1 + M3P2 - M2P3 - M1P4 - M8P5 - M7P6 + M6P7 + M5P8 = 0, 

M3P1 - M4P2 - u\Vz + M2P4 - M7P5 + MsPe + M5P7 - MePs = 0, (102) 

M2P1 - M1P2 + M4P3 - M3P4 + MePS - M5P6 + M8P7 - M7P8 = 0. 

If we wish to solve this problem from the point of view of the Hamilton- Jacobi equation 
we use the relation 

^v^ 2 A^^ fe p_ 

2,=i ''■ ^' v^s? + si + si + si + si 



z^A:=i ^fe I ^^ e=i ^ e=i ^9 ^ 

The corresponding Hamilton- Jacobi equations are 

I h^r - i, 12 /\ 2 2 -E = o, (103) 

2 j=l dSj R ^sl + g2 ^ g2 ^ ^2 ^ g2 
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' hl!-f-ii^,+E)j:u]+^{E '^ 



2D^frCdu 



D' 



Ua 



D^ 



(104) 



subject to the constraints 



dS dS dS dS dS dS dS dS 

Ui-R hM3^ U2-^ Ui- Ms- U-j- VUq- V U^-^— = 0, 

OUi OU2 OU3 OU4 OUf^ OUq OU-j OUs 

dS dS dS dS dS dS dS dS 

Uz-R M4^ Ui- hM2^ U^- hMs^ h ^5^ ^6^— = 0, 

OUi OU2 OU-i OUi UUr^ OUq OUj OUs 

dS dS dS dS dS dS dS dS 

U2-R ^1^ hU4^ ^3^ hM6^ «5^ h Ms^ ^7^— = 0. 

OUi OU2 OU3 OU4 OU5 OUq OU7 OUs 

3.2.2 Quantum motion 

The Schrodinger equation for the S-dimensional quantum Coulomb problem 



-Af )vl> +{E + ^ , '' 

2 ^ ^Jsj + sl + sl + sl + sl 



^ = 



transforms to the 8-dimensional oscillator equation (see Appendix) 



1 (8w f, tu'D^ 1 



^ = (2v-t;i 



with constraints 

T,<l> = 0, 
where operator T is given by formula ( |155| ), 

u'D' = 2 (d'E - 2z/i + ^) , 

and 

in complex spherical coordinates 



6 



Considering the oscillator equation 
(see Appendix) 



^^ — e^*^ sm^ X JT 
sm X OX OX 



uj'D^-tSD^- 



-■ + — ^ 

smx sm X 



$ = 0, 



where the operator M^ has the form 



M^ 



1 d . 3 d L^ 

sin^^a^ '''' de sin^^ 



J2 



2 6 ' 
COS"' I 



(105) 



(106) 
(107) 

(108) 

(109) 

we get 

(110) 



(111) 
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and 



J = L + f, J^ = L'^ + f^ + 2L-f. (112) 

As before, we make the complex transformation (|39| ) and also complexify parameter fi by 



putting k = ifi. We make the separation ansatz [ITS 



$ = R{i^)Zi9)G{a, P, 7; an, Ph, 7h) 



(113) 



where G is an eigenfunction of operators L^, T^ and J^ with eigenvalues L{L+1), T(T+1), 
J{J + 1), respectively. Correspondingly the wave function Z{6) is the eigenfunction of 
operator M^ with eigenvalue A(A + 3). Because there is L ■ T interaction the eigenvalue 
equation 



f2 



J G{a, /3, 7; an, Ph, 1h) = J{.J + l)G{a, /3, 7; a^, /3^, 7/^) 



(114) 



can not be separated in variables (a, /3, 7; a^if , /3//, '-yu) but we can apply the rules for the 
addition of angular momenta L and T and, following |T^ express G as a Clebsch-Gordan 
expansion 

GZ-^Tt = E {J,M\L,m'- T,t') P^,^,(a,/3,7) AS(«/.,/3^,7//) (US) 

M=m'+t' 

where {JM\Lm;Tt) are the Clebsch-Gordan coefficients. Note that the functions Gj^^.rp^ 
satisfy the normalization condition 



I ail I dilHG]^ rp^Gy^i rpn, 
Jn Juh 



27r2 



2n' 



2L + 1 \2T + 1 



SjJ' ^LL' ^TT' ^MM" ^mm' ^W ■ (US) 



If we substitute ansatz ( |113D into the Schrodinger equation ( |11U| ), then after separation 
of variables we obtain the differential equations 



1 d . n ^dZ 



A(A + 3) 



2L(L + 1) _ 2J(J + 1) 
1 — cos 9 1 + cos 9 



Z = 0, 



;ii7) 



1 d 7 dR 

^^r^ 7^ sm 9-— + 
sm V dv dv 



2n4 



(2D^e + .^-D'). ''^'^'^ ^° 



with real parameters 



£ 



2^ ^D^y 



sin d 



uj^D^ = 2{D^E-2k + 



COS^-i? 
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:ii8) 



(119) 



Consider equation ([1171) . Taking the new function by v{9) = {sin 9) 2 Z (9) we obtain 
the Poschl- Teller equation. Then the solution Z{9) = Z(^{9) orthonormalized by the 
condition 



Z'^''{9)Zir{9) sin:' 9 d9 = 6, 



\X' 



(120) 



22 



has the form 



zi^{o) 



(2A + 3)(A + J + L + 2)!(A - L - J)\ 



\ 22^+2^+2(A - L + J + 1)!(A - J + L + 1)! 
X (1 - cos^)'' (1 + cos^)^ Pi2^+^'2^+i)(cos 6), ne = 0, 1, 2, ... 



ng 



(121) 



where A is quantized as \ — L — J = rig. 

Let us now turn to the quasiradial equation (|118|) . Setting w{{}) = {sm'd)~2R(^{})^ we 
can rewrite this equation in the Poschl- Teller form 



~dtP 



4 sm v cos^ w 



w = 0. 



(122) 



Solving this equation we have following expression for quasiradial functions R{'d) 



RnAO) = (sin^)2^(cos^)'^+5, 

X 2i^i(-n^r, n^ + z^ + 2A + 4; 2A + 4; sin^^), 71^ = 0,1,2,... (123) 

with energy levels given by 
1 



£ 



2D' 



■[{n + l)(n + 8) + (2z/ - l)(n + 4)], n = 0, 1, 2... 



(124) 



where f = [uj'^D^ + 7)^5 and principal quantum number 



4/ ' 



n = 2{nr + A) = 2{nr + ng + L + J). 

Thus, the full wave function $ is the simultaneous eigenfunction of the Hamiltonian 
and commuting operators M^, J^, L^, T^, J3, L3 and T3. The explicit form of this function 
satisfying the normalization condition (see Appendix) 



iD^ 



327r2 Js^c 



^JLT ^JLTo V- 2 



dv{u) 



Uc 



IS 



$ 



JLT 

rirXMrnt 



'(2L + l)(2T+l) 
CnA^) — 1^^ RuA^) Zi\9) Gi'l^^a, P, 7; an, Ph, 1h) (125) 



where Rnr\{'&) is given by formula ( [123| ) and 
4 



CnA^) 



iiy{iy + 2X + 2nr + 4)r(2A + z/ + n^ + 4)(n^ + 2A + 3)! 



(126) 



(2A + 3)! \ Di%2(i _ e2*-<^)(A + n, + 2)(n,)!r(z/ + n, + 1) 
Let us construct now the five-dimensional Coulomb system. The constraints tell us 

f'^^u)=T{T+l)(^{u) = 0. (127) 
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and therefore the oscillator eigenstates span the states with T = and L = J. For L = J 
the Jacobi polynomial in ( p.21| ) is proportional to the Gegenbauer polynomial [^| 



pt\l''''^'\cose) 



(2/ + l)!(2L + A + 2)! ^a-2lIcos( 



:i28) 



and we obtain 



zi\e) 



Zix{9) 



22^+1 r(2L + ^) 



(2A + 3)(A-2L)! 

7r(A + 2L + 2)! 



\2L ^2L+ 



sin^)'^ Cn2(cos0). (129) 



Then from properties of Clebsch-Gordan coefficients (JM|Lm;00) = SjiSMm' and using 
l^ooi^H, Ph,1h) = 1 we see that the expansion ( |115| ) yields 



'-^Lm;00> 



a, p, 7; an, Ph, 1h) = T>^^m' («> P, l) ^JL him'- 



(130) 



Thus, the function $ now depends only on variables ("(9, 6', q;,/5,7). Observing that A 
ng + 2L = 0, 1,2, ...n, introducing the new principal quantum number A^ = (n^ + A) 
I = 0, 1, 2.. and setting k = ifi, we easily get /^from the oscillator energy spectrum (^2_ 
the reduced system energy levels 



E 



N 



N{N + 4) 
ZR2 



/i- 



2(A^ + 2)^ 



(131) 



Noting that u = ia — {N + 2) and taking into account the formulas (|123|) and ( |125| )-( p!30|) , 
we finally have the solution of the Schrodinger equation ( |105| ) as 



^ 



L\ 



Ax,0;a,f3,^) = dI el^^<!>^'{x,e-a,f3,^) 



N^^{a)R^^x{d)ZUO) 



I2L + 1 



where Zl\{6) is given by (|129|) and 



Rnr\{x) 



Kxi^) 



in x)^ e-'""^^-^-'"^ 2Fi{-N + \\ + 2 + 



sm 



2A+2 



e 2 



[(iV + 2)2 + a2](iV + A + 3)! 



(2A + 3)!\ 2R^Ti{N + 2){N -\)\ 



2^2 ^ln.'("' 


A 7) 


(132) 


- ia; 2A + 4; 1 - 


-e^^'^). 


(133) 


r(A + 2 + 2a) . 




(134) 



Thus, we have constructed the wave function and energy spectrum for the five-dimensional 
Coulomb problem. In the contraction limit R -^ oo for finite A^ we get the formula for 



the discrete energy spectrum of the five-dimensional Coulomb problem [^] 

,2 



lim En{R) 

i?— >oo 



/^ 



{N + 2) 



AT = 0,1, 
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Taking the limit R ^ oo and using asymptotic formulas as in ([56| ) we get from (|132| )-( p^ ) 



^lim ^^^^^^,(x,^;a,A7) = RNx{r)ZLx{9) J^^V^^^,{a,P,j) (135) 



iFi(-iV + A;2A + 4; ^^"^ 



nth 






R^rJr'] - 




(A^ + A + 3)! f 2fir y^ e ^^+2 


j^Nxyi ) 


(AT -A)! \N + 2) (2A + 3)! 



A^ + 2' 



which coincides with the five- dimensional Coulomb wave function obtained in paper |37 



4 Summary and Discussion 

In this paper we have constructed a series of mappings 820 -^82, S^c -^^3 and Sgc —^85, 
that are generalize those well known from the Euclidean space Levi-Civita, Kustaanheimo- 
Steifel and Hurwitz transformations. We have shown, that as in case of flat space, these 
transformations permit one to establish the correspondence between the Kepler-Coulomb 
and oscillator problems in classical and quantum mechanics for the respective dimensions. 
We have seen that using these generalized transformations ([T^) , (|6^) and (|9^) we can 
completely solve the quantum Coulomb system on the two-, three- and five- dimensional 
sphere, including eigenfunctions with correct normalization constant and energy spectrum. 

For the solution of the quantum Coulomb problem, first we transformed the Schrodinger 
equation to the equation with oscillator potential on the complex sphere. Then, via com- 
plexification of the Coulomb coupling constant fi [fi = Ze^) and the quasiradial variable 
X this problem was translated to the oscillator system on the real sphere and solved. 

It is interesting to note that the complexification of constant Ze^ /R and the quasiradial 
variable were first used by Barut, Inomata and Junker p^ in the path integral approach 



to the Coulomb system on the three-dimensional sphere and hyperboloid, and further were 
applied to two- and three- dimensional superintegrable systems on spaces with constant 
curvature [^ 0. The substitution used in p4 



e^>^ = -coth/5, /5g(-oo,oo) (136) 

is correct as an analytic continuation to the region < Rex < tt and — cxd < \m < and 
translates the Coulomb quasiradial equation with variable x to the modified Poschl- Teller 
equation with variable 13. It is possible to show that there exists a connection between 
( |136| ) and generalized Levi-Civita transformations on constant curvature spaces. Indeed, 



for instance, along with the mapping S2c^S2 we can determine a mapping H2C— ^82, i.e. 
from the two-dimensional complex hyperboloid to the real sphere: 

2 I 2 I 2 / 2 2 2\2 

Si + S2 + S3 = (M3 -u^- U2) . 



This transformation has the form 



si = t\/ul-ul-ul- — 



M? — Uo 



2U3 
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S2 = «V^3-^*l-^2 , (137) 



2 2 I Ui -\- U2 



S3 = V"3 -Ml -M2 • M3 



2m,. 



and translates the Schrodinger equation for the Coulomb problem on the sphere to the 
oscillator problem on the complex hyperboloid. Then the substitution ( |13(j| ) transforms 
the oscillator problem lirom the complex to the real hyperbolid, a solution well known 
from papers [^, |2^ . 



The method described in this paper can be applied not just to ( |TT| ) but to many 
Coulomb-like potentials. In particular the generalized two-dimensional Kepler- Coulomb 
problem may be transformed to the Rosokhatius system on the two-dimensional sphere 

As we have seen, in spite of the similarity of transformations (^ and (|T^ on the 
sphere and Euclidean space there exist essential differences. Equations (p!2D, ([6^ ) and 
(P^) determine the transformations between complex and real spheres or in ambient spaces 
a mapping C2p+i — ^ Rp+2 for p = 1,2,4. Evidently these facts are closely connected to 



Hurwitz theorem [^, according to which the nonbijective bilinear transformations satisfy 
the identity 

sj + sl + .... + s} = {uj + ul + ... + ulY (138) 

only for four pair of dimensions: (/, n) = (1, 1), (2, 2), (3, 4) and (5, 8), which corresponds 
to a mapping i?2p — > Rp+i for p = 1,2,4 respectively. 

For transformations between real spaces of constant curvature the situation is more 
complicated, and more interesting. For example, the two-dimensional transformation on 
the hyperboloid is 

2 



si = Jui ±ui±ui 



Ui -U2 
2U3 ' 



.2 = V«i±«?±ui-^ (139) 

^3 



S3 = \/ui±ui±ui- iu3± 2^ 



and 



si - si - si = {ul±ul±uiy. (140) 

Thus, the upper and lower hemispheres of the real sphere or the upper and lower sheets 
of the two-sheet hyperboloid in w-space map to the upper and lower sheets, respectively, 
of the two-sheet hyperboloid in s-space. 
The next example is the transformation 



2 



■^1 = yuf + u^-u^ 



2 2 ^1 ^2 



2u 



3 



S2 = Juj + ul-ul^^^, (141) 



S3 = \/ui + u^-u 
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2 1 „,2 „,2 I „, ""l + ^^2 



2u. 



and 



2 I 2 _ 2 
Si -\- $2 S3 



(«1 + «2 



?-7,?^2 



n 



3; 



(142) 



Here the one-sheet hyperboloid in u-space maps to the one-sheet hyperboloid in s-space. 
From transformations (|139| ) and (|141|) (using methods as in §2) it is easy to show that in 
the contraction hmit D ^ 00 this transformation goes to the real Levi-Civita transforma- 
tion (up to the translation Uj -^ V^Ui) (^.This shows that the method of this article can 
be adapted to treat a Kepler- Coulomb system on the two- and one sheet hyperboloids. 

Finally, note that in this article we do not discuss two important questions. First is the 
correspondence between integrals of motion for Kepler- Coulomb and oscillator systems. 
Second is the connection between separable systems of coordinates (not only spherical) 
under mappings (0), (0) and (0). This investigation will be carried out elsewhere. 



5 Appendix 

We present some differential aspects of the generalized Levi-Civita, KS and Hurwitz 
transformations. These calculations we are related to those in |T^, 113 for flat space. 



5.1 Transformation 820 —^82 

The Laplace-Beltrami operator on the u - sphere in complex spherical coordinates ([1 
is 

1 



A?) 



D2 

2i 



{Uidu2 - U2du,) + {U^du2 - U2du-,) + {u-idu2 - M2C?«3)' 



1 d 



— sin VG \ sin v — ~h 

D'^ 1 sin X dx dx sin^ x ^^^ 



d^ 



(143) 



while the usual Laplace-Beltrami operator on the s - sphere in spherical coordinates (x, '■p) 
has the form 



Al2) 



i?2 

1 



{sids^ - 52^.1)^ + {s^ds2 - S2dszY + {s^ds., - S2<9, 



^S-i, 



1 d 



d 
■smxT^ 



1 92 



B? \ sin x^X ^ dx ' sin^ x '^V'^ 
The two Laplacians are connected through 



A(2) 



A(2). 



The volume elements in u and s - spaces are 



dv{u) 



-e^^dxdif, dv{s) = R sinxdxd^p 



and 



-dv{s) 



ui + ui , , , 

5 — dv[U). 

ui 



(144) 



(145) 



(146) 



(147) 
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We have (the variable y? runs the from to An) 



/ ....av[s) = — — / .... 2 — dv[u). 

J So. Z J Son 



Us 



'148) 



5.2 Transformation S4C -^ S3 

The Laplace-Beltrami operator on the u - sphere in (x, q;,/5,7) coordinates is 



Al^) = §sinxe-^>^ 



-e ^ sm Xtt- 



sin xdx 



dx sin X 



(149) 



where 



/ d d cosa d ' 

Li = ? cosacot/iJ- hsmaT— -. — -7— 

\^ oa op sm p (77 ^ 

/ ^ (9 9 sina 9 ^ 

L2 = qsmacot/y- costtT-- ^stt 

\ oa op sm p (77 ^ 

^3 = -«^-, 
(7Q; 



(150) 



and 



(92 



+ cot P—- + 



d_ 1 / d^ 
df3 ^ sin^ (3\d^ 



^d d (92 

- 2 cos /3 7— 7— + 



.(9/32 Qp sin /3v 972 d^ da da^ 

while the usual Laplace-Beltrami operator on the s - sphere in (x, /5, o;) coordinates is 



(151) 



i?2 



d 



d 



1 



(92 



sin2 X (9x <9x sin2 x \ (9/32 

The two Laplace-Beltrami operators are connected by 



cot p-— 



1 92 



8(3 sin2/5(9a2 



2i 

D2 



Ar-* = ^ sin xe 2 ^ 



^^Al3) + fl-.cotx)+ .^ ■',. 
4 / sm X sm p 



^ir hr ~2cos/37— 

(97 V (97 ack , 



e2 



^/3 



and the operator acting on functions of variables (x, (3, a) is 



A?) 



-^5 



Ut 



ul + ul + ul + ul 



•aw 



2 + --^ ^-^^ 



2\ n 



D2 « /}4^ 4 ^g 

The volume elements on S4C and S3 are given by 

dv{u) = — — e2*'^sinxsin/3rfxrf/3dQ;(i7, dv{s) = R^ sin^ x^^^PdxdP da, 



u-\il52) 
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where 






dv{u) = ^^ dv{s) d-f. 



(153) 



Integration over 7 G [0, Ait] gives 



.dv(s) 



S3 



27r 754c 



ul 



dv{u). 



5.3 Transformation Sgc^Ss 

The Laplace-Beltrami operator on the u - sphere in {xy^'^'^y Pyly^^H, Ph,1h) coordi- 
nates is 



1 



D^ I sm X <^X (^X sm x 



1 9.3,9 
sm V- 



sin -(9 c^-i? d-d 



4(L^ + 2L-fsin^f + f^sin^f; 
sin^ "(9 



(154) 



where operator L is given by ( |150D and T is 

Ti = ? cos a J/ cot /^/h 7^ \-smaHT7^ ^-5-^ — 

V dan o(3h sm (3h O'Jh 



( n 9 ^ smau d 

T2 = I \ sm an cot Pht:; cos an- 



da 



H 



dPn sin PHd-fH, 



(155) 



T. 



d 



da 



H 



The Laplace-Beltrami operator on the five dimensional sphere in (x, ■&; a, (3, 7) coordinates 
is 



A(5) 



i?2 



1 d 



d 



sm X— 1 — : 



1 d . . d 4L2 
sm v-— — 



sin X ^X "^x sin x l sin ''^ f^"*^ 9-t} sin "(9 



(156) 



The Laplace-Beltrami operators are related by 



A(8) 



2i 



^smxe 



2 " 



1 2L-T + T^ 
4 'V sin^ X cos^ | 



D^A(5) + (^-6icotx 



e2 



2 -*i- 



and the operator acting on a function of variables (x, ^', a, P, 7) is 



A (5) _ jn ^9 



:^^^^-^(^^+f4s^^ 



U( 



(157) 
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The volume elements on Sgc and S5 have the form 

dv{u) = -8D^e'^'''sm^xsm^0dxdedndnH, 
dv{s) = R^sin^x^^'^^(^dxd6dQ, 



where 



We have 



dQ = - sin pdadpd'-y. (158) 



1 ^ 16i 

-J2u^^dv{u) = —dv{s)dnH (159) 

""9 i=l ^ 



and integration over the variables {an, PhjIh) gives the formula 



^ f \ iD^ f A 2 dvju) 
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